Collisionless modes of a trapped Bose gas 






o 






X3 

o 



> 

in 

o 

O 

o\ 

On 
-I— > 



I 

O 

o 



X 



M.J. Bijlsma and H.T.C. Stoof 

Institute for Theoretical Physics, University of Utrecht, 

Princetonplein 5, 3584 CC Utrecht, The Netherlands 

We calculate the excitation frequencies of the low-lying modes of a trapped Bose-condensed gas 
at nonzero temperatures. We include in our calculation the dynamics of the noncondensed cloud, 
and find agreement with experimental results if we assume that in the experiment both the in-phase 
and out-of-phase monopole modes are excited simultaneously. In order to explore whether this is 
indeed the correct explanation, we also calculate how strongly the modes couple to a perturbation 
of the external trapping potential. 



I. INTRODUCTION 

Since the experimental realization of Bose-Einstein 
condensation in magnetically trapped *^Rb, ^Li and ^"^Na 
|0-g[ gases, various properties of these Bose-condensed 
systems have been studied both experimentally and the- 
oretically. A challenging problem that has attracted 
much attention, but is still only partially resolved, is 
the temperature dependence of the excitation frequen- 
cies of the collective modes of a harmonically confined 
Bose-condensed gas. In principle, it's solution requires 
solving the coupled dynamics of the condensed and non- 
condensed parts of a highly inhomogeneous interacting 
gas @. 

In general, the dynamics of an interacting many-body 
system can be classified as being either in the hydro- 
dynamic or in the collisionless limit, and the nature of 
the collective excitations of the system changes as one 
goes from one limit to the other pj. In the hydrody- 
namic regime the mean free path of the quasiparticles 
is small relative to the wavelength of the collective ex- 
citations. Therefore one can assume the system to be 
in local equilibrium. In the case of a Bose-condensed 
gas, this important assumption leads to a description of 
the coupled dynamics of the condensate and the non- 
condensate clouds by the well-known Landau two-fluid 
equations. The solution to these equations have been 
studied by several authors P-p[ , and also experimentally 
one is trying to probe this regime |£|. In contrast, in 
the collisionless regime there is no local equilibrium be- 
cause the mean free path of the quasiparticles is much 
larger than the wavelength of the collective modes. Up 
to now, most experiments with Bose-condensed gases are 
performed in this regime. What is of prime interest at 
present is that these experiments, performed in an axially 
symmetric trapping potential, indicate a strong temper- 
ature dependence of the excitation frequencies for tem- 
peratures relatively close to the critical temperature Tc 
|lO| . The first attempts to explain this temperature de- 
pendence theoretically were unsatisfactory for a number 
of reasons. 

Far below the critical temperature, the collisionless 
modes of a trapped Bose gas are accurately described by 
the Gross-Pitaevskii equation. This equation is a non- 
linear Schrodinger equation that describes the time evo- 



lution of the condensate wave function at temperatures 
where the noncondensate density is negligible. It is non- 
linear because it includes the mean-field interaction of 
the condensate with itself. Theoretical calculations solv- 
ing this equation [|ll|-[19|, are in good agreement with 
measurements of the low -lying excitation frequencies in 
this temperature regime pO| , pi| . 

However, at higher temperatures the noncondensate 
fraction becomes large. Therefore, one has to include 
into the nonlinear Schrodinger equation also the effect 
of the mean-field interaction of the thermal cloud with 
the condensate. The low-lying excitation frequencies pre- 
dicted by the resulting nonlinear Schrodinger equation 
have been found numerically and show almost no tem- 
perature dependence ||22| , p3| . The lack of temperature 
dependence can partly be cured by including in the ef- 
fective two-particle interaction the many-body effects of 
the surrounding gas on the collisions. This causes the 
interaction to become strongly temperature dependent 
p3 and the frequencies of the low-lying modes then also 
depend on temperature [ p5[ . The frequency of the mode 
with azimuthal quantum number m — 2 found in this way 
agrees with experiments quite well. This is, however, not 
the case for the m = mode. 

The latter is related to a fundamental problem with 
the nonlinear Schrodinger equation used, which describes 
the dynamics of the condensate in the presence of a static 
noncondensed cloud. As a result it violates the general- 
ized Kohn theorem, which states that for a harmonically 
confined many-particle system there are exactly three 
center-of-mass modes with excitation frequencies equal to 
the three trapping frequencies. Clearly, this is caused by 
the fact that we also have to describe the time evolution 
of the thermal cloud. Therefore, we proposed to describe 
the collective excitations in the collisionless regime by a 
nonlinear Schrodinger equation for the condensate wave 
function, which is coupled to a collisionless Boltzmann or 
Vlasov-Landau equation, describing the dynamics of the 
noncondensed cloud pO,M. The resulting theory con- 
tains the Kohn modes exactly. 

Even above the critical temperature, solving the colli- 
sionless Boltzmann equation for the collective modes of 
the gas is not straightforward. As an illustration, we 
note that this equation also contains the hydrodynamic 
modes as a special class of solutions. However, we are in 



particular interested in temperatures relatively close to 
the critical temperature, where the discrepancy between 
theory and experiment is largest. Therefore, we want to 
formulate a manageable theory that has the correct be- 
havior both near zero temperature and near the critical 
temperature. With this objective in mind, we can treat 
the thermal cloud in the Hartree-Fock approximation, 
because near the critical temperature the mean-field in- 
teraction is small compared to the average kinetic energy 
of the noncondensed atoms. Moreover, near zero temper- 
ature the presence of the thermal cloud is unimportant. 

Finally, to solve the resulting set of coupled par- 
tial differential equations that describe the dynamics of 
the trapped Bose-condensed gas in the time-dependent 
Hartree-Fock approximation, we apply a dynamical scal- 
ing of the ideal gas results for both the condensate wave 
function and the Wigner distribution describing the non- 
condensed cloud. 

The paper is organized as follows. In Sec. |l|we present 
the general theoretical framework that describes coUi- 
sionless dynamics of a trapped Bose-condensed gas. In 
Sec. Ill we show how the excitation frequencies of the 
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low-lying modes can be determined by means of a dy- 
namical scaling ansatz. We give a physical motivation 
for our ansatz and argue in particular that it is accurate 
near the critical temperature. In Sec. ^\\ we show how 
the linear response to an external perturbation of the 
trapping frequency is calculated in our approach. This 
calculation is performed for reasons that become clear 
later on in Sec. |V| when we compare our results with 
the experiment. We end the paper in Sec. 
summary of our conclusions and an outlook. 



II. COLLISIONLESS DYNAMICS 

In the Heisenberg picture, the quantum-mechanical 
evolution of a many-body system of Bose particles in an 
external trapping potential can be described by a field 
operator ip{-K.,t) and a Hamiltonian 



H = dx 



V'nx,i)Hoi^(x,t) + 



- /dx' V-t(x,f)V;t(x',t)nx-x')^(x',t)V'(x,t) 



, (1) 



where T^(x — x') is the two particle interaction, the single 
particle Hamiltonian Tin is given by 



Wo = +Kxt(x), 



2rn 



(2) 



and Vexti^) = J2i 



j/2 denotes the external trap- 



ping potential. The field operator satisfies the equal-time 
commutation relations 



^(x,i),v3(x',i) 



= 0, 



V^(x,i),Vnx',i) =<5(x-x') 



(3) 



By definition, the time evolution of the field operators 
is determined by the Heisenberg equation of motion 



' dt 



V'(x,i),iJ 



(4) 



From this equation we can now easily derive an equa- 
tion of motion for the condensate wave function ^(x, i), 
which is defined as the expectation value (■0(x, t)). Ap- 
plying a mean-field approximation results in a nonlin- 
ear Schrodinger equation for 'I'(x, i), which includes the 
mean-field interaction with the noncondensate density 
n'{x,t). It reads 



ih 



9*(x,f) 
dt 



2m 



+ Vext{^) 



+ T^^ [2n'(x, t) + no(x, t)] } *(x, t) . (5) 

Here, the condensate density no(x, i) is defined as the 
modulus squared of the condensate wave function, i.e., 
|^(x, i)p. The factor of two difference between the con- 
densate and noncondensate mean-field interactions is un- 
derstood physically by realizing that the noncondensate 
contributes both a Hartree and a Fock term, but the con- 
densate only a Hartree term. Furthermore, the two-body 
interaction has been renormalized to a hard-core poten- 
tial T^^(5(x — x'), where the two-body scattering matrix 
rp2B _ ^^^ a/m solves the Lipmann-Schwinger equation 
for the scattering of two particles with zero momentum 
and a denotes the interatomic scattering length. In prin- 
ciple, the anomalous average has renormalized the two- 
body interaction potential to the many-body T-matrix 
j^MB ^ which also includes the effect of the surrounding 
gas on the collisions between two particles ||2^ . For sim- 
plicity, however, we will use the two-body scattering ma- 
trix instead, neglecting the effective temperature depen- 
dence of the two-body interactions. 

Next, we want to derive an equation of motion for the 
noncondensed cloud. This is achieved by considering the 
Heisenberg equation of motion for the one-particle den- 
sity matrix. 



(V>t(x,t)V3'(x',i)) (^'(x,i)^'(x',i)) 
(^'t(x,i)^'t(x',t)) (^'(x,t)Vi't(x',t)) 



(6) 



Here, the field operator describing the noncondensate is 
defined according to iip'{x,t) = i}{x,t) — vl'(x, t). 

From Eq. (m we derive a Boltzmann equation. This 
Boltzmann equation describes the time evolution of the 
quasiparticle distribution function i^(k, x, t), which is re- 
lated to the Fourier transform of the above one-particle 
density matrix by a local Bogoliubov transformation. 
More precisely, it is derived from the equations of mo- 
tion for the one-particle density matrix by performing 
a gradient expansion in the center-of-mass coordinate 



(x + x')/2. This is justified because in general the non- 
condensate density profile varies on a much larger length 
scale than the external trapping potential. The resulting 
Boltzmann equation contains two steaming terms, corre- 
sponding to the local group velocity of a quasiparticle and 
the local force on a quasiparticle. The velocity and the 
force are given by the momentum and the spatial deriva- 
tive of the energy dispersion -E(k, x, t), respectively. The 
distribution of quasiparticles can also change because of 
collisions and on the right-hand side there is an asso- 
ciated collision term. In total the Boltzmann equation 
therefore reads 



dE d 



d 



dt dhk dx 



dE d 
dx dhk 



F = 



OF 



dt 



(7) 



If we treat the quasiparticles in the Popov approxima- 
tion |2Sf| , their energy in the frame where the superfluid 
velocity Vs{x, t) is zero, is given by 



^(x,k,i) 



2u2 



2m 



+ yext(x) + 2r^^n(x,i) 



[T^^no{x,t 



n9{x,t) 

2x1/2 



mv1{x,t) 



Kk-^fs{x,t). (8) 



Here 

tion, i.e., ^(x. 



X, t) is the phase of the condensate wave func- 
t) 



y/no{x,t) exp i9{x, t) . Furthermore, 
the superfluid velocity is proportional to the gradient of 
that phase or more precisely Vsix, t) — hV6{x, t)/m. 

In the temperature region of interest it is sufficient to 
treat the quasiparticles in the Hartree-Fock approxima- 
tion. The energy in the laboratory frame is then given 
by 



2u2 



^(x,k,i) 



2m 



+ Kxt(x) + 2T"^ji(x,i) 



(9) 



The resulting nonlinear Schrodinger and collisionless 
Boltzmann equation are expected to describe the dynam- 
ics of the trapped Bose gas near the critical temperature 
Tc, because in this temperature region the mean-field in- 
teraction of the condensate is small compared to the av- 
erage kinetic energy of the noncondensed cloud, implying 
that the Hartree-Fock approximation is valid. Moreover, 
they will also be accurate near zero temperature, because 
here the noncondensate density is negligible. 



III. VARIATIONAL APPROACH 

Directly solving the set of equations that describe the 
coupled collisionless dynamics of the condensate and the 
noncondensate, i.e., the nonlinear Schrodinger equation 
given by Eq. dq) and the collisionless Boltzmann equation 
following from Eq. m) , is difficult even when treating the 



quasiparticles in the Hartree-Fock approximation. How- 
ever, it is well known that a simple scaling ansatz for the 
condensate wave function gives the correct frequencies of 
the low-lying modes at zero temperature [ [l4|jl^ ]. There- 
fore, we use a similar method to solve our set of equations 
and also assume the time dependence of the quasiparticle 
distribution function to be given by a dynamical scaling 
ansatz. We thus get 



no{x,t) = 






*. 



n 



(1) 



(x,i) = E 



moji 



(3) 

Vi X^ 



-P^XJ 



F(k,x,t) = 
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^/^ 


[V^^ 


^i 


kv 


m OLi 
h a.i 



Vi 



(2) 



(10a) 



(10b) 



(10c) 



-V, 



(2)i 



m .(2) 



The six parameters {Ai} and {ai} describe the coupled 
monopole and quadrupole oscillations of the two compo- 
nents of the gas. The six parameters {77- } and {r]\ '} 
are included to describe the in-phase and out-of-phase 
dipole, or Kohn modes. Notice that the momentum ar- 
gument of the scaling ansatz for the distribution function 
contains explicit time derivatives of the scaling parame- 
ters. Together with the quadratic form of the phase, this 
ensures that the continuity equation is to be satisfied. 
The variational parameters {77; } and {(3i} are therefore 

related to {tjI } and {A;} in a way that is specified be- 
low. Due to the mean-field interaction, the equilibrium 
values of {Ai} and {a^} are in general not equal to one. In 
order to ensure that the equilibrium profile is isotropic in 
momentum space, the factors {q} are inserted, which are 
equal to 1/af , where cXi denotes the equilibrium value of 
ai . The modes of the noncondensed cloud described by 
this ansatz, are however not constrained to be isotropic 
in momentum space. Hence they describe more compli- 
cated dynamics than hydrodynamic motion, as desired 
in the collisionless limit. 

The equations of motion for the scaling parameters of 
the condensate are most easily found from the lagrangian 
density for the nonlinear Schrodinger equation written in 
terms of the density and the phase, 

/: = ;.Mx,t)M!^-,^^[^"°(-'*)]^ 



dt 



2m 4rio(x, t) 



2m 



no(x,t)[V0(x,t)]2 



T 



2B 



K.t(x) + 2T^^n'(x, t) + ^— no(x, t) 



(11) 
no(x,t) . 



Inserting the scaling ansatz for the density and the appro- 
priate phase into this lagrangian density, and integrating 



over space yields a lagrangian for the scaling parameters. 
The equations of motion for these parameters are then 
given by the appropriate Euler-Lagrange equations. The 
equations of motion for the scaling parameters of the non- 
condensate arc found by taking moments of the collision- 
less Boltzmann equation with respect to ki,Xi, kikj, Xiki 
and XiXj. 

The resulting equations of motion for the in total 12 
variational parameters are 



••(2) , 2 (2) 

m + w, f]] 



2T^B d 
mN' aj^(2) 



{^n,{i^))nc\{^ , (12b) 



_T^2 r^^(ci>„„(x)), 1 -p. 1 
Ai + ujj Ai — — , o, TT H — ^^^ "^ 1 1 T~ 



'm'{x^)c Af 2m{x^)c Xi ^^ \j 



2T^^ d 
m{x'j)c dXi 



(*n'(^l)>cn— ' (12c) 



2 _ t;'^"'"' 2 ^ r2^($„'(x))„, c, T-rV^ 



ai +uj^ai^ 3 



TT v:2 



2t2^ a 

m{xf)nc dai 



(<I'no(^2)>ncnr ' ^^^^^ 



whereas the dynamics of the phase is determined by 

ZWi Ai 



respectively. Finally, the weighted averages with respect 
to the equilibrium condensate and noncondensate densi- 
ties are 



(/(x)). 



dx$„„(x)/(x) 



and 



(/(x))„, -y"rfx$„,(x)/(x) 



(17) 



(18) 



where $„'(x) = {2t:Y^ j d\i $F(k,x). Note that the 
equations of motion for {Ai} and {qj} are of the same 
form, apart from a factor of two difference between the 
mean-field interaction of the condensate with itself and 
the mean-field interaction of the noncondensate with it- 
self. Of course, this is the same factor of two difference 
that was mentioned above in relation to the nonlinear 
Schrodinger equation. 

In the Thomas-Fermi limit we can neglect the average 
kinetic energy of the condensate relative to the mean- 
field interactions, and the ground-state density profile of 
the condensate is approximately an inverted parabola. 
The opposite limit, where the mean- field interaction is 
much less than the average kinetic energy, results in a 
gaussian density profile. For the experiments of interest 
here, a gaussian profile is appropriate near the critical 
temperature. At low temperatures this is no longer true 
and an inverted parabola is more accurate. However, for 
our purposes it is convenient to take at all temperatures 
a gaussian ansatz for the condensate, because it is known 
to give the correct frequencies even at zero temperature 
PI. Hence, 



.„.„.„„ n(^) 



1/2 _y- rr^^2 



(19) 



LOi 



(13b) 



For notational convenience we have defined the argu- 
ments ^2 Etnd ^2 by 



6,i = {X^X^ -I- 1)1' - T]l')^ilai , 

6,i = (azXi/^i + rif' - r]f^)/Xi . 



(14a) 



(14b) 



Also, the kinetic energy of the condensate and noncon- 
densate are given by 



-.fem.c / j_. ^^ [^^<^no{^)? 



Tr^'-= dx 



and 



r^kin,nc^ /dx 



2m 4$„„(x) ' 
dk Ti^fc? 



(15) 



(27r)3 2m 



$j.(x,k), (16) 



Since the effect of the noncondensate is most impor- 
tant near Tc, we take the dynamical scaling ansatz for 
the quasiparticle distribution to be a Bose function. Near 
the critical temperature and neglecting mean-field inter- 
actions, this indeed gives the correct mode frequencies, 
i.e., {2uJi}. We thus take 



*F(k,x) = 



N'UAf^hio,) ^(n^e ^mujfxV 



-N\ 



-E- 



(20) 



C(3) " \ 2m ' ^ 2 

where N^e) is the usual Bose-distribution function 

N{e)=[eP^-l\-' . (21) 

Here, (3 = l/ksT and C(3) « 1.202. The normalization 
factors are such that the total number of condensed and 
noncondensed atoms is equal to Nq and N' , respectively. 
We take the temperature dependence of the total num- 
ber of condensed atoms to be given by the noninteracting 



result, ]Vo = [1 - {T/TcY]N, where N ^ N^^ + N' \& the 
total number of particles. 

The equations of motion for the scalin g p arameters 
that result from inserting Eqs. (|l9| ) and ( P0|) into Eq. 
( p^ ) are listed in appendix ^ If we introduce a vector 
notation 



n(x, t) -^ n(x) + 5n(x)e*' 



(26) 



Ul 



U2 = 



n 



A 

a 

,(1) 



u) 
a; 



(22a) 



(22b) 



(22c) 



where w^ = SijUJi, these can be schematically written as 



iii + J7^ -Ul = v(ui,U2) , 



U2 + Jl • U2 = w(ui,U2) 



(23a) 



(23b) 



with V and w nonlinear vector functions of their argu- 
ments. To find the excitation frequencies of the modes, 
we have to linearize these equations around the equi- 
librium value (ui,U2) and subsequently determine the 
eigenvectors and eigenvalues of the linearized problem. 
From the explicit expressions in appendix H, it is easily 
seen that the linearized equations for Sui and (5u2 decou- 
ple. The resulting excitation frequencies for the in-phase 
and out-of-phase monopole, quadrupole and dipole, or 
Kohn modes are presented and discussed in section M 
below. 



IV. LINEAR RESPONSE 

In order to compare our theoretical results with exper- 
iment, it turns out to be important to understand what 
is measured experimentally when the external trapping 
potential is perturbed to excite the collective modes of 
the Bose-condensed gas. It is possible that a periodic 
modulation of the trapping frequency will excite more 
than one mode when it is not exactly on resonance. The 
question is then which mode is most likely to be seen ex- 
perimentally. To answer this question we have to study 
the linearized response of the gas to a periodic perturba- 
tion of the trapping frequencies 



LUi —f uji + Sujie^'^ , 
leading to 

14.t(x) ^ l^e.t(x) + (514,t(x)e^'^* 
and 



(24) 



(25) 



A quantity that characterizes the response to such a per- 
turbation is the time averaged work done by the pertur- 
bation lEol, 



W^^Jd^SVe^t{^)5n*{x.), 



(27) 



where the asterix denotes the complex conjugate. An 
explicit expression for the time-averaged work is found 
by linearizing the density profiles around equilibrium, by 
putting Xi — Xi + ^Aic"^* and ai = cti + Saie'^'^*. If we 
insert the resulting expression for Sn ~ Shq + 5n' into 
Eq. (|^), we get 



W = }^muji5uji [Xi5X*{x1)c + ai5a*{x'l)nc\ 



(28) 



To calculate the work done by the perturbation, we thus 
need to know the response of the scaling parameters to 
a perturbation of the external potential. 



After linearizing the equations of motion Eq. (23a), 
with U2 = 0, to first order in SuJi and (5ui, the resulting 
equation of motion for the fluctuation 5\ii , reads 



- uj'^Sui + rj2 • (5ui ^ [Vn^v] • (5ui - 2Vt ■ 5Vt ■ Ul . (29) 

The partial derivative of v with respect to Ui is to be eval- 
uated in the equilibrium point. These linearized equa- 
tions of motion are easily solved by 



5\l^ 



E' 



.(") 



(30) 



,(») 



where u^ ' denote the normalized eigenvectors of the ho- 
mogeneous part of Eq. (p9) and a„ = 2\i]^ ■ n • JH • Ui . 
The time averaged work done can now be expressed as 

bn 






where the residue 6„ is given by 



6„ = ^ muiSujian 



K>^(^c 



..(")/ 



(31) 



(32) 



In reality, the eigenmodes are damped. Therefore, in 
a theory that includes damping, the poles Un have an 
imaginary part and W is always finite with a maximum 
at w = oJn- Hence, a resonance would occur if the sys- 
tem were driven with that frequency. A measure for the 
strength with which this mode is excited, is the residue 
hn- Also, the time-averaged work would acquire an imag- 
inary part, which would determine the power absorption. 
We have calculated the residue b for two particular modes 
as a function of temperature and the results are presented 
in the next section. 



V. RESULTS 



In this section we present the resuhs of our calculations 
and if possible compare them with experimental data. In 
principle our approach determines 12 modes of the gas. 
In an axially symmetric situation, they correspond to the 
in-phase and out-of-phase versions of three Kohn modes, 
two monopole modes and one quadrupole mode. These 
in-phase and out-of-phase modes are the coUisionless ana- 
logue of the hydrodynamic first and second sound modes 
P-|[. We have calculated the excitation frequencies of 
these modes for the parameters of the experiments with 
**^Rb (To), 23Na i] and ^Li |). Furthermore, for the ex- 
periments with Rb, we have calculated the residue b 
mentioned in the previous section for the in-phase and 
out-of-phase monopole modes. 

From the results presented in Fig. n^, it is clear that 
the temperature dependence of the out-of-phase m — 2 
mode is in reasonable agreement with that of the m — 2 
mode of the * Rb experiments. Furthermore, the ex- 
perimental data for the m = mode goes to the correct 
non-interacting limit near Tc^ where it coincides with our 
theoretical curve for the in-phase m = 2 mode. In addi- 
tion, between 0.7 and 0.6 Tc the experimental data drops 
to the zero-temperature limit (10/3)^/^ Wr [|ll|, where it 
coincides with our theoretical curve for the out-of-phase 
m — 2 mode. It is important to note that in obtaining 
these results, we have included the effect of evaporative 
cooling by fitting the total number of particles to the 
experimental results of Jin et al. pO| . Qualitatively, it 
appears that the strong temperature dependence found 
experimentally might be due to the fact that one simulta- 
neously excites both the in-phase and out-of-phase m = Q 
modes [Q . To explore this possibility we have calculated 
the residue h for these two modes. As shown in Fig. 0, 
there is a clear crossover at about 0.5 Tc, after which the 
value for the residue of the in-phase m = Q mode shoots 
up. Therefore, the experimental data might actually be 
due to the excitation of two modes, and as a function of 
temperature one crosses over from exciting mainly one 
mode to exciting mainly the other. Furthermore, Fig. S 
shows that the out-of-phase m = mode cannot be seen 
experimentally below 0.2 Tc, because here the mode is 
very difficult to excite. Finally, we note that it might be 
possible to observe the two additional modes present in 
our calculation experimentally. Whether this is possible 
depends on the overlap of these modes with the applied 
perturbation, and on the damping of the modes, which 
we have neglected . 

As mentioned above, our theory reproduces the cor- 
rect excitation frequencies near the critical temperature 
and near zero temperature. In principle, sufficiently far 
below the critical temperature there is a 'dimple' in the 
noncondensate density profile due to the presence of a 
condensate. This 'dimple' is not taken into account in 
our dynamical scaling ansatz. However, it is incorrect to 
include this 'dimple' by simply modifying the scaling pro- 



file, because we are using a one parameter scaling ansatz, 
Consider for example an out-of-phase mode of the con- 
densate and the noncondensate. 
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FIG. 1. The in-phase and out-of-phase m = 0,2 modes as 
a function of T/Tc for the experimental conditions of Jin et 
al. [10]. The relative phase of the density profiles of the 
condensed and noncondensed atoms is denoted by <j). Also in- 
cluded are the experimental results for the m = (triangles) 
and m — 2 (circles) modes found in these experiments. 
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FIG. 2. The residue b for the modes m = 0, = (solid 
line) and m — 0,(j} = n (dotted line) as a function of T/Tc for 
the experimental conditions of Jin et al. [10]. 

Because the 'dimple' is caused by the mean-field inter- 
action with the condensate, it also has to move out-of- 
phase with the exterior part of the noncondensed cloud. 
This implies that there are two length scales that deter- 
mine the dynamics of the noncondensate. Hence, mod- 
ifying the scaling profile ^p, for example by taking the 
exact equilibrium solution for the noncondensate density, 
can actually give worse results. Moreover, neglecting the 
'dimple' in the density profile of the noncondensed cloud 



is certainly correct near the critical temperature. A qual- 
itative reason why the theory still seems to agree well 
with experiment, is that a 'dimple' in the noncondensate 
density on the one hand increases the effective mass of 
the condensate, which reduces the excitation frequency, 
but on the other hand lowers the mean-field interaction, 
which increases the excitation frequency. Apparently, 
these two effects almost cancel each other. 

In Fig. ^ we show the calculated in-phase and out- 
of-phase monopole and Kohn modes for the experimen- 
tal parameters of Stamper-Kurn et al. |9[|. The Kohn 
modes are exactly present in our theory, which can be 

seen by rewriting Eq. (|l2| ) in terms of rji + tj^ 
(1) „(2) 



V, 



and 
77j- ' . The out-of-phase Kohn modes have the qual- 



itatively correct feature found in the experiments, that 
the frequency of the mode is shifted downward with re- 
spect to the trapping frequency. Experimentally the shift 
is about 5%. Instead, we find a reduction of about 10%. 
The fact that the out-of-phase dipole mode becomes un- 
stable at T/Tc « 0.05 can be understood by realizing 
that our equilibrium profile is always centered around 
the origin. Because we do not include the 'dimple' in 
the noncondensate density profile this implies that at a 
certain temperature, it is energetically favorable for the 
noncondensed cloud to shift it's center outward. Because 
the instability occurs only at very low temperatures, this 
artifact of our ansatz is not of importance for our pur- 
poses. 
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FIG. 3. The in-phase and out-of-phase 7n = (solid lines), 
and Kohn (dashed lines) modes as a function of T jTc for the 
experimental conditions of Stamper-Kurn et al. [9], with a 
total number of particles A'^ = 40 x 10^ 

To be able to compare our data with results found in 
the Popov calculations of Zaremba et.al. 122], we have 
also looked at the case of an isotropic trap. Interestingly, 
our results, which are shown in Fig. 0, are quite similar 
to theirs. Moreover, also when analyzing the tempera- 
ture dependence of the mode frequencies by means of a 
temperature dependent effective interaction, the results 



found in the isotropic case are very similar to these re- 
sults |32|. This suggests that the interesting temperature 
dependencies are related to the anisotropy of the external 
trapping potential. 
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modes as a function of T/Tc for an isotropic trap 
ujx ~ (.^y = i^z ~ ^0 = 200Hz and with a total number of 
atoms N = 2000. 
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FIG. 5. The in-phase and out-of-phase I = 0,2 modes 
(solid lines) and Kohn (dashed lines) modes at constant tem- 
perature T — SOOnK, as a function of the number of atoms 
in the condensate, for the experimental conditions of Bradley 
et al. [2]. The condensate collapses for Nq = 1447. 



Finally, as shown in Fig. ||, we have also calculated 
the mode frequencies for the experimental conditions of 
Bradley et al. |g] that apply to the case of a negative 
scattering length. It should be noted that our approach 
is particularly suited for a discussion of this case, be- 
cause the mean-field interaction of the condensate is at 
most comparable to the energy splitting in the trap, due 
to the intrinsic instability of the condensate to collapse 



p3[ . As expected, the maximum number of condensate 
particles is slightly shifted downwards by the mean-field 
interaction of the noncondensate with the condensate. 
At a temperature of 300 nK , we find a shift of about 
3%, whereas more accurate calculations find a shift of 
about 5% [Q. In addition, the dynamical treatment of 
the noncondensed cloud seems to have little effect on the 
mode frequencies. This indicates that it is a good ap- 
proximation to treat the collapse purely in terms of the 
condensate even at nonzero temperatures |35 . 



VI. CONCLUSION AND OUTLOOK 

If we interpret the temperature dependence of the 
m — Q mode found in the JILA experiments as due to the 
excitation of two modes instead of one, there seems to be 
reasonable agreement between our theoretical predictions 
and the experimental results. We should mention, how- 
ever, that our results are sensitive to the explicit form of 
the ansatz. Indeed, work done by the authors in collabo- 
ration with E. Zaremba, shows that using for the ansatz 
a numerical solution for the condensate and nonconden- 
sate equilibium density profiles in the Popov approxima- 
tion, leads to qualitatively different results. However, as 
mentioned above, we believe this to be due to the one- 
parameter nature of the ansatz used. It is a matter of 
future investigation to resolve this problem, possibly by 
including more variational parameters into the density 
profile. Alternatively, we can perform a RPA- calculation 
of the mode frequencies Q , which is essentially equiva- 
lent to finding the modes of the collisionless Boltzmann 



equation coupled to the nonlinear Schrodinger equation 
exactly [g . Finally, we can in principle also calculate the 
damping of these modes within our approach, by includ- 
ing the collision terms in the Boltzmann equation [p7[ . 
Work to implement these ideas is in progress. 
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APPENDIX A: EQUATIONS OF MOTION 



We list here the equations of motion for the varia- 
tional parameters {Ai}, {a^}, {rjl '} and {77^^ '} result- 
ing from inserting the explicit form of the condensate 
wave function and the noncondensate distribution func- 
tion Eqs. ^ and ^ into Eqs. (|l|). 



The average 
harmonic oscillator length is defined by I = ^h/mu!, 
where lj = {u}i0J2(jJ3)^'^ ■ The thermal wavelength is given 
by Ath = j2TTn'^/mkBT and C(4) w 1.082. Note that 
in equilibrium 77^ = 77, = and that the equations 



earized around this equilibrium 



for {Xi,ai} and {7?j , t?,-^''} decouple when they are lin- 
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